Casimir Energy and Entropy in the Sphere-Sphere Geometry 
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We calculate the Casimir energy and entropy for two spheres described by the perfect metal model, plasma 
model, and Drade model in the large separation limit. We obtain nonmonotonic behavior of the Helmholtz free 
energy with separation and temperature for the perfect metal and plasma models, leading to parameter ranges 
with negative entropy, and also nonmonotonic behavior of the entropy with temperature and the separation 
between the spheres. This nonmonotonic behavior has not been found for Drude model. The appearance of this 
anomalous behavior of the entropy is discussed as well as its thermodynamic consequences. 



I. INTRODUCTION 



H. MULTIS C ATTERING FORMALISM OF THE CASIMIR 
ENERGY 
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In 1948, Casimir predicted the attraction between per- 
fect metal parallel plates (11 and between neutral polarizable 
atoms IH due to quantum fluctuations of the electromagnetic 
field. Some years later, Schwinger extended this formalism to 
dielectric plates at finite temperature Jll]. Recently, a multi- 
scattering formalism of the Casimir effect for the electromag- 
netic field has been presented J3,|5t] (see also |0] and 10]). 

The Casimir effect has some peculiarities. In particular, it 
is a non-pairwise interaction; the Casimir thermal force (the 
thermal part of the Casimir energy) between two isolating 
bodies is not necessarily monotonic in their separation, as seen 
in the sphere-plate and cylinder-plate cases H. In addition, 
for some geometries, intervals of negative entropy appear, as 
in the case of two parallel plates described by the Drude model 
|0] or, as recently shown, in the interaction between a Drude 
model plate and sphere ifioll and in the interaction between a 
perfect metal plate and sphere ifTTI . 

In this article we study the Casimir effect between two 
spheres in the large separation approximation using differ- 
ent models for the electric permeability: (1) the perfect metal 
model, (2) the plasma model, and (3) the Drude model. As a 
result, we find negative entropies in certain ranges of temper- 
ature and separation between the spheres for the perfect metal 
model, and for low penetration length for the plasma model. 
In addition, we find nonmonotonic behavior of the entropy 
with the separation while the force is attractive for all separa- 
tions, making it appear as though the natural evolution of the 
system tends to increase the entropy in certain ranges of tem- 
perature and separation. For long plasma length and for the 
Drude model we do not find this anomalous behavior of the 
entropy. We also discuss the thermodynamical meaning and 
consequences of negative entropies in Casimir effect. 

The remainder of the article is arranged as follows: In Sect. 
Hfl we describe the multiscattering model used herein to ob- 
tain the Casimir energies and entropies for the two spheres. In 
Sect, [ill] we obtain the large separation limit of the Casimir 
energy between perfect metal spheres, as well as the entropy 
and force. We also obtain entropies at smaller separations nu- 
merically. We study the plasma model system in Sect. HVIand 
the Drude model in Sect. [V] We discuss the thermodynamic 
consequences of these results in Sect. [VI] Finally, we discuss 
the results obtained in the Conclusions. 



To calculate the Casimir energy, entropy, and forces be- 
tween two spheres, we employ the multiscattering formalism 
for the electromagnetic field fl, qJ]. This formalism relates 
the Casimir interaction between objects with the scattering of 
the field from each object. The Casimir contribution to the 
Helmholtz free energy at any temperature T is given by 



E = k B Tj2 Iog|I-N(«n)| 



(1) 



n=0 



where n„ = are the Matsubara frequencies and At = 
2 Jj, c rp is the thermal wavelength. The prime indicates that 
the zero Matsubara frequency contribution has height of 1/2. 
All the information regarding the system is described by the 
N matrix. For a system of two objects, this matrix is N = 
T1U12T2U21. Ti is the T scattering matrix of the i& object, 
which accounts for all the geometrical information and elec- 
tromagnetic properties of the object. is the translation ma- 
trix of electromagnetic waves from object i to object j, which 
accounts for all information regarding the relative positions 
between the objects of the system. 

For a sphere of radius R with electric and magnetic 
permeabilities e and /1, the T matrix is diagonal in 
(imP,£'m' P') space, with elements given by T^ f , TO , = 
Tf m , where the Tf m are defined as H 
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£K e+ i(x), and 



where i t {x) = yjl^ifi), k e {x) = 
n = ^/Jj£- Expressions for the U Q/ 3 matrices can be found in 
H and El. 



III. PERFECT METAL SPHERES 

In the perfect metal limit, we apply e — > 00 for any /x to 
Eqs. (O and ©. In this case, the T matrix elements are inde- 
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pendent of /i, taking the well-known universal form 

"2 K t+ i(KR)' 
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A. Casimir energy in the large separation limit 

To obtain the large separation limit of the Casimir energy, 
we need the dominant part of the T matrix in this limit. We 
define the adimensional frequency q by k = q/d, where d is 
the distance between the centre of the spheres. The main con- 
tribution in the large separation limit comes from the lowest- 
order expansion term of the T matrix elements in l id. As it is 
known that, at small n, the T elements scale as n 2(+1 0], the 
dominant contribution comes from the dipolar polarizabilities 
part of the T matrix, taking the form 
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By the use of the universal relationship between determinants 
and traces, log \ A\ = Trlog(yl), and applying a Taylor expan- 
sion in terms of i of Eq. (fl3, we obtain 



E = -k B Tj2 Tr (N(A^ 1 n)) 



(7) 



Using the translation matrices in a spherical vector multipole 
basis and the large separation approximation of the T ma- 
trix, the trace of the N matrix is obtained by straightforward 
calculus. Here we denote with a sub-index T the results valid 
for all temperatures, with a sub-index the results in the 
quantum limit (T — > 0), and with a sub-index cl the results 
in the classical limit (Ti 0), which is equivalent to the high 
T limit. Carrying out the sum over Matsubara frequency, we 
obtain the Casimir contribution to the Helmholtz free energy 
as 
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(2 (15 - 29z 2 + 99z 4 ) cosh(z) + 15cosh(5z) 
+ (-45 + 58z 2 + 18z 4 ) cosh(3z) 
+24z (6z 2 - 5 + (5 + 3z 2 ) cosh(2z)) sinh(z)) , 



(8) 



where z = d/Xr- We also define the adimensional Casimir 
energy as E a d{z) = j^eEt- From this result, the quantum 
(T —> 0) and classical (h — > 0) limits with their first correc- 
tions are easily obtained as 
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where z = d/Xr- Note that in Eq. (O, the first correction 
to zero temperature case is p roportional to T 6 , contrary to 
the plate-sphere case JH ll ill and to the cylinder-plate case 
JH], where a result proportional to T 4 were obtained for both 
systems. The result presented here is new for perfect metal 
spheres, but expected, because it has been obtained also for 
compact objects in the diluted limit lfl3i [Till . As correc- 
tions have the same negative sign of the Casimir energy in 
both limits, they describe an increase of the magnitude of the 
Casimir energy in high and low temperature limits. It is no 
longer the case if we study the energy for all temperatures. 
In fact, for some ranges of separation and temperature, ther- 
mal photons tend to reduce the Casimir energy between the 
spheres, as shown in Fig. [TJ This effect is captured in the 
contribution of the next term to the quantum limit, this term 
is AE Q = + ^dknT^R 6 (ksZy an d tends to reduce the 
magnitude of the energy. Fig. [TJis an indicator of the appear- 
ance of negative entropy in this system because of the negative 
slope of the energy curve. 
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Figure 1. (Color online) Casimir energy between perfect metal 
spheres as a function of compared with the energy at zero tem- 
perature. The dotted curve is the quantum limit, the dashed curve 
is the classical limit, and the solid curve is the asymptotic finite- 
temperature Casimir energy. Note that these curves are independent 
of the radius of the spheres but are only valid in the large separation 
limit. 



It is clear in Fig. Q]that for some temperatures and distances, 
the Casimir energy between spheres is lower than in the zero 
temperature case. But it is less evident the validity Eq. ((9). In 
fact, there is a tiny increase of E/ En at low until reaching 

a local maximun at ^ w 1.0388 of E/E as 1 + 1(T 4 . It is 
not visible in Fig. Q]because the difference os scales. The local 
minimum is reached at ss 3.21733, where E/En ~ 1 — 
3.46 x 10 -2 . This minimum is reached, in the important cases 
of room temperature (T = 300K) and at the boiling point of 
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N 2 (T = UK) at d T = 
respectively. 



300 



3.89/ito and g?t=77 = 15.21/im 



B. Casimir entropy in the large separation limit 

In the canonical ensemble, the entropy is defined as S = 
—8tE, where we remain that we denote the Helmholtz free 
energy by E. In the large separation limit, the Helmholtz 
free energy depends on the adimensional variable z = = 

2tt Jt T , so we can write the entropy as 



he 



S 



dz dE 
dT~dz~ 



-2-k 



dk B dE 
He dz 



(11) 



and define the adimensional entropy as S a d(z) = ^ 
—d z E a d(z). From this result, the quantum (T — > 0) and clas- 
sical (h — > 0) limits are easily obtained as 



So = 
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Because of the limit at low temperature of the entropy, we 
know that the entropy is positive for low T (not seen in Fig. 
|2]because it is small compared with S/S c u but it can be ob- 
served in Fig. [3}, so there are three points where 5 = 0, 
including the origin. These two new points where the en- 
tropy nulls correspond to the local maximum and minimum 
observed in the Casimir energy in Fig. Q] 

Negative entropy of the Casimir effect has already been ob- 
tained between Drude parallel plates in Jgl and in Il5ll . and 
more recently between a perfect metal plate and sphere in lUIll 
and between a Drude sphere and plate in lfl6ll and ifioll . This 
is the first time, to the best of our knowledge, that negative en- 
tropy appears between spheres because of the Casimir effect. 

These results are only valid when the separation between 
the spheres is large compared with their radius, regardless of 
the radius of each one. Therefore, it is possible that the inter- 
val of negative entropy would disappear when the separation 
between the spheres reduces, regardless of the temperature of 
the system. A numerical exploration of entropies at smaller 
separation has been performed to verify this issue. 



C. Numerical study at smaller separations 



15k B R 6 
4d 6 
R 6 
2d 6 



(13) 



^— ( 15 + 30z + 27 z 2 + Uz 3 + 9z 4 - 18z 5 ) e~ 2z , 



where z = j^. So, the entropy is a growing function with 
temperature in both limits, but this is not the case for all tem- 
peratures (the next low temperature expansion term of Eq. 

(EJi is AS Q = - ls ™?f k B dR 6 (^) 7 , which already indi- 
cates that the entropy could change its growing behavior with 
the temperature), as we can observe in Fig. [2] where a region 
of negative entropy and another region of negative slope of the 
entropy are observed. 
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Figure 2. (Color online) Casimir entropy between perfect metal 
spheres as a function of compared with the classical limit. The 
dashed curve is the classical limit, and the solid curve is the asymp- 
totic finite-temperature Casimir entropy. Note the regions of negative 
entropy and negative slope of the entropy with the parameter -A- 



As noted in the previous subsection, asymptotic results are 
no longer valid when the separation between the spheres be- 
comes comparable to their radius. For this reason, a numerical 
study of entropy was performed for these cases. We computed 
Eq. (Q~|i numerically for all temperatures from T = until 
reaching the classical limit for fixed ratio between the radius 
R and separation, r = 4. 

For spheres, the T matrices are diagonal but infinite ma- 
trices, so a cutoff in (£, £') space is needed to obtain finite- 
dimensional matrices. In addition, another cutoff in Matsub- 
ara frequency is needed to obtain a finite series. 

The proposed method is an asymptotic approximation, at 
small separations more and more modes are needed to obtain 
convergent results. This means that there exists a minimum 
separation between the spheres below which we are not able 
to take into account enough multipoles to obtain good results. 
We use multipoles up to £ < 15, which means that we are 
restricted to r max = -A- < 0.45, when contact is reached at 
- R — _R 

2R 



0.5, where the energy diverges. 



In the small separation limit, the proximity force approx- 
imation (PFA) is known to be a good approximation to the 
Casimir energy It is also known that perfect metal plates 
do not experience negative entropies, so we do not expect to 
observe negative entropies between spheres in the small sepa- 
ration limit. 

In Fig. [3] the entropy of the system of two perfect metal 
spheres is plotted as a function of z — for constant r = ^ . 
The large separation result is shown too. We choose a log- 
log representation of the absolute value of the entropy divided 
by its corresponding classical limit. Therefore, zeros are ob- 
served as log divergences, and we can also observe the cases 
of negative entropy. Starting in the large separation regime, 
we observe an interval of negative entropy. As we increase r 
(reducing the separation between the spheres), the region of 
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negative entropy tends to reduce until it disappears between 
r = 0.40 and r = 0.41. Power-law decay of the entropy (TTZt 
at low temperatures is observed as a linear decay of the curve 
at low z, and constant behavior in the high-temperature limit 
( fT3] l is also reached in the computation. 
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temperature force is plotted as a function of -y-. Here, for 
constant temperature, we observe nonmonotonic behavior of 
the force with the adimensional parameter . 

Nonmonotonic behavior of the adimensional Casimir force 
implies nonmonotonic force behavior with temperature, but 
not with separation, because of the extra dependence of the 
force on the separation in Eq. (TPfl i. In fact, it is easy to verify 
that the force behaves monotonically with separation, and the 
nonmonotonicity of the entropy with separation implies non- 
monotonic behavior of the force with temperature, because 



dF 
df 



d 2 E 
dTdd 



dS_ 



(15) 



so the appearance of negative slopes of the entropy with sepa- 
ration implies nonmonotonicity of the Casimir force with tem- 
perature, despite the attractive force for all separations and 
temperatures. 



Figure 3. (Color online) Log-log plot of the absolute value of the 
entropy divided by its classical limit for perfect metal spheres at con- 
stant r as a function of Starting from the asymptotic solid blue 
curve (r —5- 0), we increase r in steps of 0.1 up to r = 0.4 (magenta 
curve for r — 0.1, yellow for r — 0.2, red for r = 0.3 and r = 0.4 
is the dotted curve). Dashed curve is the case r = 0.41, where the 
interval of negative entropy has disappeared. The curve for r — 0.45 
(solid green curve) is also shown 



The power law decay of the entropy reduces from S oc T 5 
for asymptotic results (r — > 0) to around S oc T 3 for the 
closest studied case (r = 0.45). 

D. Casimir force in the large separation limit 

In the previous subsection we demonstrated that intervals of 
negative entropy appear due to the Casimir effect between per- 
fect metal spheres. In these cases, for any given temperature, 
a minimum of entropy exists for a given, nonzero separation. 
Naively, this would imply a violation of the third law of ther- 
modynamics if the Casimir force is not zero at the minimum 
of the entropy. However, this is not the case, as will be ex- 
plained in Sect. [Vl] In addition to that, we find that the force 
is always attractive, independent of the increase or decrease of 
the entropy, which also would naively imply a violation of the 
second law, because the system can be enforced to perform 
a process in which the entropy is reduced instead increased. 
However, this is not the case, as will be also explained in Sect. 
IVII However, the appearance of negative entropy does have an 
effect on the force. The asymptotic Casimir force F = —ddE 
can be written in terms of the adimensional Casimir energy as 

F ad {z) = ^4fi F = 7E ad (z) - zd z E ad (z), (14) 

he Ry 

where z = In Fig. |4j the adimensional asymptotic force 
between the perfect metal spheres compared with the zero- 
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Figure 4. (Color online) Asymptotic force divided by its zero- 
temperature limit as a function of . The dotted curve is the zero- 
temperature result. The dashed curve is the classical limit, and the 
solid curve is the result at finite temperatures. The nonmonotonic be- 
havior of the force with temperature results from the negative slope 
of the solid curve at constant separation 



As observed in Fig. |4] the force between the spheres is 
always attractive, but it is not monotonic with temperature; 
asymptotically, for any given temperature, there exists a range 
of separations for which the force decreases with temperature. 
This is the first time, to the best of our knowledge, that non- 
monotonic behavior of the Casimir force with temperature has 
been described between compact objects. Nonmonotonicities 
of the force between a plate and cylinder and between a plate 
and sphere were already obtained in |j8|], but in that case the 
nonmonotonicity already appears for the scalar field. Non- 
monotonicity does not appear between spheres for the scalar 
field; this is a characteristic effect of the electromagnetic field, 
because cross-polarization terms of the Casimir energy are es- 
sential for this nonmonotonicity to appear. 
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IV. PLASMA MODEL 

In this section, we assume that the electric susceptibility of 
both spheres is described by the plasma model, i.e., 



e(icn) 



1 



(2tt) 2 

XpK? 



(16) 



and fi = 1. To obtain the large separation limit of the Casimir 
energy, we need the dominant part of the T matrix in this limit. 
The main contribution comes from the dipolar polarizabilities 
part of T matrix, taking the form ifTTl [l6ll 



r 



MM 
lm,lm' 



RXW 
I2ir 2 d 3 



(3 + y 2 -3ycoth(y)) 



R lm,lm' 



qRY 
d 



(17) 



(18) 



with y = 2tt-^. The coefficients of the magnetic polariz- 
ability depend on the plasma frequency, tending to the perfect 
metal result as Xp and to zero in the transparent limit 

X P -)• oo HI- 



A. Casimir energy in the large separation limit 

Once we have the asymptotic T matrix, the Casimir energy 
can be obtained by a straightforward but long calculation. It is 
possible to obtain analytical results for finite temperatures, but 
they are too long to show here. The zero-temperature Casimir 
energy is 



E 



hcR 6 



l6Ttd 7 y 4 

(207 + 222y 2 + U3y 4 - 414j/coth(y) 
+207y 2 coth(y) 2 - 222y 3 coth(y)) , 

and the high-temperature limit is given by 



(19) 



3i? 6 
~~d? 



1 + ^4 (3+y 2 -3ycoth(y)) : 



(20) 



metal limit, obtaining negative entropy and nonmonotonic be- 
havior of the entropy with separation and temperature as be- 
fore. However, when Xp 3> R, the spheres are transparent to 
the magnetic field and the problem reduces to a scalar prob- 
lem. In this case, the entropy is always positive. In Fig. [5] the 
adimensional asymptotic entropy, compared with its classical 
limit, is plotted as a function of for these two limits. In 

Fig. |6j we present the points ^pf) corresponding to zeros 
of the entropy (continuous curve), its temperature derivative 
(dashed curve), or its separation derivative (dotted curve). As 
discussed in the case of perfect metal plates, the anomalous 
behavior of the entropy found here would naively imply a vi- 
olation of second and third laws of thermodynamics, because 
the Casimir force is always attractive, irrespective of the slope 
of the entropy. However, this is not the case, as we will dis- 
cuss in Sect. [VT] As noted above, in the perfect metal limit 
we obtain a region of negative entropy, while for low A p this 
region disappears. The point where negative entropies appear 
depends on the magnetic susceptibility [i. For fi = 1, negative 
entropies appear around ApK 2R. 
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Figure 5. (Color online) Entropy between plasma spheres as a func- 
tion of ^ divided by the classical limit for the two limit cases. The 
solid curve is the entropy in the perfect metal limit Xr <C R, where 
a region of negative entropy and nonmonotonicities of entropy with 
separation and temperature are present; the dashed curve is the en- 
tropy in the transparent limit \r 3> R, where the usual monotonic 
behavior of entropy is shown 



with y = 2n The perfect metal limit is reached for Xp — > 
0, as expected. When the plasma wavelength Xp ^> R, only 
the electric sector of the T matrix contributes, so in this case 
the energy in the zero- and high-temperature limits is given by 

23hcR e 3i? 6 



V. DRUDE MODEL 

In this section we assume that the electric susceptibility of 
both spheres is described by the Drude model, i.e., 



B. Casimir entropy in the large separation limit 



e(icK) = 1 + 



4tt 2 



(22) 



The entropy of the system is obtained in the same way as 
for the perfect metal case, but for the plasma model we have 
two different regimes. When Xp <C R, we reach the perfect 



and /j, = 1. To obtain the large separation limit of the Casimir 
energy, we need the dominant part of the T matrix in this limit. 
The main contribution comes from the dipolar polarizabilities 
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Figure 6. (5^7 , ) map of zeros of the entropy (dotted curve), the 
temperature derivative of the entropy (dashed curve), and the separa- 
tion derivative of the entropy (dotted curve) 



part of the T matrix, taking the form |[TTl[l6ll 

t mm _ Jtn R<J ( qR\ 4 rEE 2 (qR\ 3 

*-l m ,l m >- 45 - y d ), hmArn'-^ [jj ) ■ W 

Now, the response of the material changes dramatically from 
the case of the plasma model. In principle, one would expect 
to obtain the T matrix for the plasma model as a — s- 00, but it 
is obvious that this is not the case. Not only do we not recover 
the plasma model in this limit, but also the T matrix diverges. 
The reason is that the Plasma model is some kind of singular 
limit of Drude model, with a qualilative different behaviour. 
Now the dominant contribution is given by the electric sector 
of the T matrix, whereas the magnetic part does not contribute 
for asymptotic energies. 

A. Casimir energy and entropy in the large separation limit 



model, leading to an effective scalar Dirichlet problem. The 
entropy is obtained as in the perfect metal case. As only the 
electric polarization contributes, we expect the usual behavior 
of the entropy as a mono tonic function of as shown in Fig. 
[7] The zero- and high-temperature Casimir entropy limits are 
given by 

_ 7Q4k B R 6 fnk B T\ 5 11776k B dR 6 (-Kk B T\ 7 

~ 315d \ he J 1575 \ he ) 1 

(26) 

R 6 

S c i=3k B — . (27) 

In this case, the correction term to the quantum limit tends to 
reduce the entropy, but it is not strong enough to change the 
sign of it, as seen in Fig. [7] 

s/s cl 
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Figure 7. (Color online) Entropy as a function of the adimensional 
parameter divided by its classical limit for the large separation 
limit of the system of Drude spheres. Entropy is positive and mono- 
tonic in temperature and separation 



Once we have the asymptotic T matrix, obtaining the 
Casimir energy is a straightforward calculation. Carrying out 
the sum over Matsubara frequency, we obtain the Casimir con- 
tribution to the Helmholtz free energy as 

2k B TR 6 e 5z 
d e {e 2z - if X 
((6 - 10z 2 + 22z 4 ) cosh(z) + (-36s + 12z 3 ) sinh(z) 
+ (l2z + 4z 3 ) sinh(3z) + 3 cosh(5z) 



En 



+ (-9 + 10z 2 + 2z 4 ) cosh(3z)) , 



(24) 



where z = d/Xx- The zero- and high-temperature Casimir 
energy limits are given by 



E 



23hcR 6 
47rd 7 



k R T 



3i? 6 



(25) 



These asymptotic Casimir energies are universal because they 
do not depend on any material property. The lack of mag- 
netic polarizability contributions is characteristic of the Drude 



VI. THERMODYNAMICAL CONSEQUENCES 

In this section we discuss the thermodynamical conse- 
quences of the obtained results. In this article we have ob- 
tained and compared the large separation limit of the Casimir 
energy and entropy for two spheres using three different elec- 
tric susceptibility models. For perfect metal spheres, at any 
nonzero fixed temperature, we observe an interval of separa- 
tions for which entropy is negative, while at zero temperature 
the entropy is always zero. 

In addition, the Casimir force is attractive for all separations 
and temperatures. So, we could naively think that we have 
possible violations of the third and second laws of thermody- 
namics, due to the existence of processes where the entropy 
of the system tends to increase and to the negative entropy 
intervals at finite temperature and distances, respectively. 

According to the Krein formula fl7il . we know that the 
Helmholtz free energy of the electromagnetic field has three 
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independent additive contributions (see Eq. (l35l l of the Ap- 
pendix): one is from the thermal bath, being proportional to 
the volume of the space lfl8ll . Another is the sum of contri- 
butions of objects immersed in the bath considered as isolated 
objects, each contribution being also a function of the volume 
and surface of each object 111811 . The third is of geometrical 
nature, which we could call the Casimir part of the Helmholtz 
free energy and that we actually calculate in Eq. ([TJ. This 
term depends on the electromagnetic nature of the objects in 
the system, but it has a geometric nature, because it depends 
on the relative separations and orientations between the ob- 
jects, being zero iff the relative separations between all the 
considered objects become infinite. 

Considering the whole system, the nonmonotonicity of the 
Casimir entropy with temperature is compensated by the con- 
tribution of the vacuum, because one scales with the global 
volume lfl8ll (Eq. (l36l l) and the other with the separation be- 
tween the spheres (Eq. ([T3V). This is not the case for the 
nonmonotonicity with separation. At constant temperature, 
the entropies of the thermal bath and of each object remain 
constant, but the Casimir force is always attractive while the 
entropy can increase or decrease with separation. Therefore, 
the potential violation of the second law of thermodynamics 
still requires an explanation. 

The Krein formula also states that internal sources of en- 
tropy inside the spheres cannot appear in order to compensate 
these regions of anomalous behavior of the entropy, unless 
these internal sources are independent of the separation be- 
tween the spheres. 

The second law of thermodynamics states that global en- 
tropy must increase for any process, but only in closed sys- 
tems. As we are working in the canonical ensemble, we 
are implicitly assuming that there exists an external reservoir 
which keeps our system at constant temperature, so the sys- 
tem is not isolated and entropy can increase or decrease with- 
out violation of the second law. In the canonical ensemble, 
the condition equivalent to the second law is that the global 
Helmholtz free energy must decrease for any process, and 
this is true for the studied system. Therefore, the appearance 
of nonmonotonic entropy behavior in the canonical ensemble 
just implies nonmonotonic behavior of the force with temper- 
ature, as seen in Eq. (TT3T >. 



VII. CONCLUSIONS 

In this article we have obtained and compared the large 
separation limit of the Casimir energy and entropy for two 
spheres using three different electric susceptibility models. 
For perfect metal spheres, at any nonzero fixed temperature, 
we observe an interval of separations for which entropy is neg- 
ative, while at zero temperature the entropy is always zero. 

We can trace the origin of this anomalous effect to the 
functional form of the Casimir energy in this large separation 
regime. It consists of four components, two of which are of at- 
tractive nature, equivalent to the sum of two scalar problems, 
one for the electric polarization and the other for the magnetic 
polarization. The other two components come from the cross- 



coupling between the electric and magnetic polarizations of 
the two spheres, and tend to reduce the Casimir energy be- 
tween the spheres (Eq. 12 of fll an d Eq. 6 of lfl3lo . These 
cross-coupling polarization terms are responsible for the im- 
possibility of factorization of the electromagnetic Casimir en- 
ergy into two equivalent scalar problems in general, and ap- 
pear because the translation matrix U Q( g is not diagonal in po- 
larization space. In this article we show that these cross terms 
not only reduce the Casimir energy, but also for some sepa- 
rations and temperatures, their contribution to the entropy is 
greater than the contributions of direct coupling between the 
electric and magnetic polarizabilities, resulting in an interval 
of negative entropy and nonmonotonic behavior of the entropy 
with separation and temperature. 

In addition, the Casimir force is attractive for all separations 
and temperatures. So, we could naively think that we have 
possible violations of the third and second laws of thermody- 
namics, due to the existence of processes where the entropy 
of the system tends to increase and to the negative entropy 
intervals at finite temperature and distances, respectively. 

In Sect. I VII having into account the complete thermody- 
namical system, that the system is described by the canonical 
ensemble, and with the help of the Krein formula and Weyl 
formula, we have demonstrated that there are not violations of 
second and third laws respectively. Therefore, the appearance 
of nonmonotonic entropy behavior in the canonical ensemble 
just implies nonmonotonic behavior of the force with temper- 
ature, as seen in Eq. J15) . 

The interval of negative entropy for the spheres appears be- 
cause of the cross-coupling between the polarizations of the 
electromagnetic field, which leads us to conclude that it is a 
characteristic of the electromagnetic field and does not have 
an analog in the Casimir effect due to scalar fields. 

Applying the PFA to this problem, we obtain that, for per- 
fect metal spheres near contact, the entropy is always positive, 
so we performed a numerical study of the entropy between 
perfect metal spheres of equal radius at smaller separations. 
The results showed that there exists a minimum separation be- 
tween the spheres for which the negative entropies disappear. 
In addition, the region of negative slope of the entropy with 
separation disappears for another smaller separation. 

We have also obtained the energy and entropy for the 
plasma model, for which similar results are obtained for 
plasma penetration length Ap < 2R, while for Ap > 2R, 
the negative entropies disappear, and for another greater Ap, 
the nonmonotonic behavior of the entropy with separation also 
disappears. 

In the transparent limit (Ap ^> R), the magnetic polariza- 
tion does not contribute to the asymptotic Casimir effect, so 
the problem reduces to a scalar field problem and we find that 
entropy is a positive monotonic function. 

When we study the asymptotic Casimir effect with Drude 
model spheres, the system is qualitatively different. The 
magnetic polarization does not contribute to the asymptotic 
Casimir effect because it depends on d~ 4 instead of the d~ 3 
dependence of the electric polarization term. Then also for 
Drude spheres the asymptotic limit reduces to a scalar field 
problem, resulting in a positive monotonic entropy. 



8 



Nonmonotonicities of the Casimir force between objects 
are not unusual. As the Casimir effect is not pairwise additive, 
the interaction between two objects is affected by the presence 
of a third, leading to nonmonotonicities of the Casimir force 
between cylinders |[l9i l20ll . or between spheres 1I2D1 in the 
presence of a plate. The nonmonotonicity presented in this 
article has a different nature, as it does not come from inter- 
actions with a third object but rather from correlations with 
temperature, as seen in Eq. (TT5V 

In addition, in the assumptions of the multiscattering for- 
malism it is implicitly assumed that the objects are stationary, 
at fixed positions in space, so all results shown here are valid 
for quasistatic processes. Therefore, we must carefully con- 
sider when this quasistatic assumption does not apply, because 
the stationary system could abandon equilibrium, requiring 
more careful study 1I22I l23ll . 
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APPENDIX: KREIN FORMULA FROM 
MULTISCATTERING FORMALISM 

In this appendix, we will derive the Krein formula of the 
electromagnetic field fnll from the multiscattering formula 

a. 

Krein formula factorizes the density of states of the electro- 
magnetic field in three parts lfl7ll . Here we present a deriva- 
tion of such formula from the Multiscattering formalism of 
the Casimir effect flU]. 

The partition function (Z = J^^Lo °f trie m 
the presence of N general dielectrics of arbitrary geometry is 
obtained, before regularization of the Casimir energy, as Q5|] 



(28) 



\Sn\\MV 



where \S n \ 1 = | A + k„\ 1 i s me n— Matsubara frequency 
contribution to the partition function of the thermal bath, and 
M is the next N x N block matrix 



M a p = S^T- 1 + (S al3 - l)V aP . 



(29) 



If we multiply and divide Eq. d28l ) by |A^oo| (being Moo the 
same Ai matrix as above, but with each object at an infinite 
distance from each other), the Helmholtz free energy of the 
whole system can be writen as = log(Z)) 



n=0 



log|S„|+log|7W c 



log 



\M\ 

\Moo 



where log \ S n \ is the contribution of the thermal bath in ab- 
sense of inmersed objects, which leads to the Planck formula 
of the blackbody spectrum, log|A^oo| = Sa=i 1°§ l^al * s 
the contribution to the Helmholtz free energy of each dielec- 
tric object immersed in the bath considered as an isolated ob- 
ject, and log fppr^yl = l°g |I — N| is the Casimir part of the 
Helmholtz free energy, which depends on the electromagnetic 
nature of the objects in the system, but it has a geometric na- 
ture, because it depends on the relative separations and orien- 
tations between the objects, being zero iff the relative separa- 
tions between all the considered objects become infinite. This 
term is the Casimir energy that we actually calculate in Eq. 
0- 

The energy of a quantum system is related with the Density 
of States p(uj) by IU7tl 



dco— ojpioj). 



(31) 



Using Eq. d30b . the same energy in the zero temperature case 
can be writen as 



he f°° 
E= — / dn 
2Wo 



log|5„|-£log|T a |+Iog|I-N| 



a=l 



(32) 



Because the dispersion relation of the massless EM field d 
ck, and integrating by parts, we obtain 



r 00 hi 
E = I duj-Lu—x 

2 7T 



(33) 







-d u log 1^1 + J2 5 - lo § l T «l - 9 - lo S I 1 - 



Then, because Eqs. (T3TT > and d33l ), we can relate the Density 
of States with scattering properties of the system as 



p(ui) = - 



-d u \og\S u \ + Y2 log |T a I - log |I - 



(30) 



(34) 

As a conclussion, we can factorize the Density of States as the 
sum of three independent terms 

n 

/?M = PoM + 5^PaH +pc{u), (35) 

where po{oj) = —\du log \Sq\ is the Density of States of the 
free EM field, p a [iS) = i<9 w log |T a | is the change of the 
vacuum Density of States of the EM field because each a— 
dielectric object considered as a isolated object in the medium, 
and pc{u) = —^du log |I — N| is the change of the Density 
of States related with the relative positions of the dielectric 
objects, which is zero when all the objects are infinitely apart 
from each other. pc( w ) is me source of the Casimir effect. 

Eq. d35l ) is the Krein formula 1 17i], derived from the Multi- 
scattering formalism of the Casimir effect for the EM field. In 
fnll . Eq. (03 at zero temperature is obtained from the Krein 
formula removing the divergent contributions. Here we have 



followed the opposite way, a derivation of the density of states 
of the a thermal bath with intrusions and it factorization in Eq. 
<|35). 



flew mm, 



V 



(36) 



Terms po(uj) and p a {oj) diverges with w because main 
contribution to Weyl formula for 3D volumes of the EM 



but pc( w ) converges if our system consists on TV compact 
objects, because in this case, N can be demonstrated to be 
a trace class operator for all cj, then the determinant is well 
defined S. 
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